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1. Introduction 

Whilst the singlet sector of QCD has been well studied using lattice QCD within the quenched 
approximation and even with Nf = 2 there have been few studies using Nf = 2 + 1. Computational 
cost is the main reason for this — the computation of correlators for singlet mesons involves the 
calculation of disconnected diagrams which are inherently noisy and so require a long Monte Carlo 
timeseries in order to be accurately measured. 

We are engaged in a study of the scalar and pseudoscalar singlet sectors of QCD with Nf = 
2+1 flavours of "Asqtad" improved staggered quarks [[!]]. The MILC ensembles have given Asqtad 
a strong track record in terms of physical results and they are one of the cheapest of the fermion 
formalisms to simulate, allowing us to build up a large number of configurations. 

2. Simulation Details 

Using the UKQCD's QCDOC [§ we have generated these large ensembles at 2 different lattice 
spacings using the one-loop tadpole-improved Liischer-Weisz gauge action [^] The generation of 
the p = 7.095 ensemble is still underway and is over halfway to completion. The lattice spacing 
has been obtained through a determination of ro/a from the static quark potential (taking ro = 
0.467 fm). From here on the /3 = 6.75 ensemble will be referred to as the "coarse" ensemble, and 



Nf J3 I? x T 


ami am s ro/a a [fm] 


Ncfg Mraj Target 


2+1 6.75 24 3 x 64 
2+1 7.095 32 3 x 64 


0.006 0.03 3.8122(74) 0.12250(24) 
0.00775 0.031 5.059(10) 0.09230(19) 


5375 32250 30000 
1911 11466 20000 



Table 1: Ensembles generated. 



the p = 7.095 ensemble will be referred to as the "fine" ensemble. 

2.1 Algorithm 

The coarse ensemble was generated using the RHMC algorithm as formulated in [Q], with a 
second order leapfrog integrator. This has the benefit of being an exact algorithm, compared to the 
inexact R algorithm [§] used by the MILC collaboration which has 0(8r 2 ) errors, and so allows 
for larger integrator timesteps. The fine ensemble also used the RHMC algorithm but exploited the 
improvements of Clark and Kennedy [§] which include mass preconditioning of the light fermion 
kernel with the strange quark, and the n f/, -root multiple pseudofermion trick. Combined these allow 
us to use the fourth order Omelyan integrator with the gauge fields and fermion fields on different 
timesteps. 

2.2 Measurement 

When measuring fermionic operators using staggered fermions we must take into account both 
the spin and taste structure of the state in question. To this end we work in the spin ® taste -basis 
notation of Kluberg-Stern. The states which we consider are the singlet scalar 1<8>1 and the singlet 
pseudoscalar 75 (8)1. In principle we could also use the 7475 <8 1 operator for the singlet pseudoscalar 
but for reasons to be discussed we use the 75 (8 1- 
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The connected diagrams are measured using standard point sources, and we use the stochastic 
volume source method to measure the disconnected diagrams. Using Gaussian noise sources 
we compute the operators for the F$ ® Fj state thus 

VrW4lI I Tn^Ar^rvM- 1 ^ , (2.1) 

s i x,X4=ty,y4=t 

where the Ar s ®r> covariantly and symmetrically shifts the source in the hypercube and applies the 
relevant Kogut-Susskind staggered phase in order to obtain the correct state. The disconnected 
contribution to the correlator is then calculated using 

Dr s ®r T (At) = (4 s ®r T (0 (t + At)) (2.2) 

with the appropriate VEV subtraction for the scalar singlet. 

For certain choices of Ar s ®r r (those that separate source and sink by an even number of links) 
we are able to apply a variance-reduction trick [f7|] which allows us to use eight or fewer noise 
sources and still obtain a stochastic error comparable to the gauge error. This is applicable to both 
the 75 g) 1 and 1 (g> 1 operators, and is why we choose 75 (gi 1 over 7475 <g> 1. However since it is only 
applicable to a subset we still use N s = 32 in order to obtain a good estimate of other operators. 

2.3 Disconnected Statistics 

As has been mentioned previously disconnected contributions carry with them the noise of 
the fermionic sea and so in order to reduce statistical noise a large number of configurations is 
required. Even more troublesome for the 75 <g> 1 is its relation to the topological charge which has a 
notoriously long autocorrelation time. 

In Figure [I] we show the reduction in relative error achieved from performing such a long run. 
This improvement should carry over directly to an improvement in the error on an estimate of the 
r]' — n mass difference since the ratio of the disconnected to the connected correlator is related to 
this. For further details on the ^-ratio and the statistics of the disconnected contribution the reader 
is referred to our methods paper [JSj] . 

3. Scalar flavour-singlet Spectroscopy 

Whilst the heavy scalar singlet mesons (XcO, Xbo) WQ we U understood in terms of their con- 
stituent quark content the light scalar meson sector is far from as well understood. The interpreta- 
tion of the light scalar singlet states — the fo resonances — has eluded many experimental studies, 
phenomenological models and even lattice studies. 

The particle data table lists five scalar singlet states with mass less than 2 GeV — the fo (600) [a] , 
/o(980), /o(1370), /o(1500) and /o(1710). The quark model predicts only two light isoscalar qq 
states (/ and /') so the picture for the light isoscalar mesons is obviously more complicated. A pop- 
ular conjecture is that the three heavier states are produced by a mixing of the two flavour singlet 
states / and /' with the light scalar glueball. 

The scalar glueball has been well studied in quenched lattice QCD and is accepted to have 



a mass in the region 1.5 — 1.7 GeV [^, [10p . However in dynamical simulations the sea quarks 
cause the glueball and ++ qq interpolating operators to couple to the same states and the glueball 
becomes a less well defined concept. 
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where e(Af ) is the standard error on the mean of D ql/ (At)) on 



Figure 1: The relative error (e(Af) = J^) . 
the light-light disconnected correlators for the 75 ® 1 (with VKVR applied, and local source and sink) for 
various subsets of the full UKQCD long coarse run, shown against the j3 = 6.76, m; = 0.01 , m s = 0.05 MILC 
coarse ensemble. 



3.1 Glueball Interpolating Operators 

In order to perform a mixing study it is necessary to design our operators so that they couple as 
strongly to the desired state as possible. In order to couple to a mainly glue state G with J PC = ++ 
we use the following operator 

tf«(p,0 = £ (oUx,t) + □(,(*, + n^fct)) (3.1) 



where the D l u are plaquettes in the kl -plane with i = 0, . . . , 3 levels of Teper blocking [11] applied 
to the gauge field, on top of 2 levels of APE smearing [[J with smearing constant c = 2.5. We have 
measured these operators at p • p = and p • p = 1 for the fine ensemble, and so far at p • p = for 
the coarse ensemble. With our set of operators we form the 2x2 correlator matrix C as 

Qj(r) = (#l(t)t?j(t + T)) - (3.2) 

for the zero-momentum operators, and 

Qj(T) = (^(t)ffj(t + T)) (3.3) 

for the p • p = 1 operators, and use the variational method to obtain a set of operators with maximum 
projection onto the ground state. Finally we use the ratio of the variational eigenvalues Ao,Aj to 



project out the remaining contamination [12], obtaining our final mass estimate 



am eff (t) - ^am eff (t - 1) 

amr\[t) = 5 . (3.4) 

1 M 
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Figure 2: rom e ff(0 vs f / r o f° r both the fine and coarse glueball operators, with both p • p = and 1 for the 
fine and p • p = for the coarse. The effective mass has been calculated using equation 3.4. 



3.2 Taste Symmetry Violation 

Rather than removing the doublers completely, the staggered formulation reduces the 16 dou- 
blers to 4 degenerate "tastes" of fermion. These are then reduced to one by taking the fourth-root 
of the fermion determinant. Unfortunately the taste symmetry is broken by interactions with the 
gauge field which leads to large mass discrepancies between hadrons of different tastes. The Asqtad 
action [|l|] has improved taste-symmetry violations (0{a^a 2 )), and is the most widely used stag- 
gered action. The physical validity of the fourth-rooting procedure has been called into question, 
but there is a large body of theoretical work confirming the validity of the resulting rooted action 



(see [13, |4!] for reviews), although not everyone is convinced [lla 



In previous analyses of the scalar non-sing let state a (980) /a { 1450) @, it was noticed 
that the mass obtained was significantly below the lowest allowed decay threshold m % + . Indeed 
its mass appeared to be more consistent with a %% state which, due to G-parity, is forbidden as a 
decay channel for the scalar non-singlet. This was explained within the framework of staggered 
chiral perturbation theory (SChPT) as being due to unphysical multi-hadron states being introduced 
by the taste-split T]-multiplet [|TJ]. The most significant of these is the Tin state which at large time 
completely dominates the ciq correlator. 

It is also the case that there are such unphysical contributions to the /o correlator, though 
fortunately there are explicit predictions from SChPT for these contributions to both the fo and ao 
correlators j jl9[ ] which in principle should allow us to include the contaminations in our fits. 
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Figure 3: Our measurements of the ++ glueball mass plotted against a 2 , with previous UKQCD analyses 
of the ++ glueball on MILC ensembles and UKQCD quenched ensembles using the Wilson gauge action. 
The horizontal lines are at the positions of the /o resonances. 



3.3 Results 

Analysis of the scalar-singlet fermionic operators is ongoing, so we present results for the glue- 
ball interpolating operators only. Our final mass estimates are obtained from the weighted average 



of the projected effective mass (equation [3.4] ) over the plateau region. For the coarse ensemble this 
is done for the momentum-zero effective mass on 4804 configurations betweeen t = 1 and t = 4, 
and for the fine ensemble for the momentum-one effective mass on 1000 configurations between 
t = 2 and t = 5 (see Figure ||). This gives masses in physical units of Mg(0 ++ ) = 1629(32) MeV 
for the coarse ensemble, and Mg(0 ++ ) = 1600(71) MeV for the fine ensemble. These are shown 
in the context of previous lattice determinations of the ++ glueball mass in Figure [| 

Currently the p • p = effective mass shows no plateau for the fine ensemble, and we have yet 
to measure the momentum-one operators for the coarse ensemble, but it is our intention that the 
weighted average be performed over both p • p = and 1 when this becomes possible. 

4. Pseudoscalar flavour-singlet Spectroscopy 

The light pseudoscalar flavour-singlet mesons — v\ and t]' — are particularly interesting due 
to the role topology plays in the mass of the latter, and the importance of disconnected diagrams in 
calculations involving them. It is also an important test of the validity of the fourth-root trick used 



with staggered fermions. It has been suggested [15] that even after the continuum and chiral limits 



have been taken the chiral behaviour of staggered fermions will not match that of QCD. Study of 
the 77' should help show whether or not this is the case. 
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The improvement in the statistical error on the disconnected operators from such high statistics 
has been shown in Figure [[} Analysis of our correlators is ongoing so we do not present any results 
here, but for a detailed account of our motivation and methodology the reader is referred to our 
exploratory study [)8j|. 
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